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Time-resolved diffraction microscopy technique has been used to observe the formation of laser-induced pe¬ 
riodic surface structures (LIPSS) from the interaction of a single femtosecond laser pulse (pump) with a 
nano-scale groove mechanically formed on a single-crystal Cu substrate. The interaction dynamics (0-1200 
ps) was captured by diffracting a time-delayed, frequency-doubled pulse from nascent LIPSS formation in¬ 
duced by the pump with an infinity-conjugate microscopy setup. The LIPSS ripples are observed to form 
sequentially outward from the groove edge, with the first one forming after 50 ps. A 1-D analytical model of 
electron heating and surface plasmon polariton (SPP) excitation induced by the interaction of incoming laser 
pulse with the groove edge qualitatively explains the time-evloution of LIPSS formation. 


Femtosecond laser-induced periodic surface structures 
(LIPSS, or ripples) have been actively studied during 
the last decade, due to applications of surface structur¬ 
ing such as laser micro/nano machining 1 , solar celld^, 
wave guideeP, super hydrophobic surfaces^, and many 
other+^. 

While the mechanism by which femtosecond pulse in¬ 
teraction leads to LIPSS formation on different materi¬ 
als is under debat^, the model for describing the period 
and orientation of LIPSS with spatial period A > A/2 
(laser wavelength A) is held in strong consensus 89 . In 
this model, the incident laser excites and interferes with a 
surface electromagnetic wave/surface plasmon polariton 
(SPP), and creates periodic intensification of the fields 
at the interface, which somehow results in groove for¬ 
mation with the same period. While SPPs can only be 
excited natively in metals, LIPSS also may be formed 
on dielectrics and semiconductors given sufficient laser- 
induced modification of the dielectric function via inter¬ 
band transition, so as to become metallic during the laser 
pulsd 9 . In addition to metallicity, SPP excitation from 
free-propagating radiation requires a surface structure or 
roughness. Femtosecond LIPSS formation is therefore 
understood essentially as a multi-pulse effect, where a 
phase-matched grating structure is developed at A only 
after of several pulses of ablation undergoing positive 
feedback!^. 

The timescale of LIPSS formation is fundamentally in¬ 
teresting for femtosecond pulses because both the laser 
and SPP (^ fs)would have long since vanished before 
the expected timescale of material movement (^ ps). 
And while dynamics of the femtosecond laser damage 
process (without LIPSS) has been studied using pump- 
probe techniquejm, the vast difference between the dam¬ 
age timescale and typical laser repetition rates (ns ~ ms) 
has limited these measurements to dynamics primarily 
of single-pulse damage. Consequently, the multi-pulse 


nature of femtosecond LIPSS makes them a challenging 
subject for a pump-probe study. So far there has only 
been one effort to temporally resolve their formation 12 
where Murphy et. al. used a double-pump pulse and 
probe, to observe some evidence of periodic structure for¬ 
mation on Si surface at +50 ps, but they were not able 
to follow the dynamics afterward and the whole interac¬ 
tion was complicated by non-uniform damage crater and 
possibly residual heating/incubation created by the first 
pump pulse. In this work, we present a method to pro¬ 
duce LIPSS with a single femtosecond pulse by way of 
a mechanically-formed groove on a copper target in air, 
and dynamically resolve LIPSS formation with Time Re¬ 
solved Diffraction Microscopy (TRDM) for the first time. 

For this experiment, single crystal Cu (MTI Corp. 
10 x 10 x 1mm, (100)) was chosen as a target instead 
of dielectric or semiconductor targets like sapphire or 
Si because the dielectric function of Cu naturally pro¬ 
vides support for SPP creation and eliminates other time- 
dependent variables such as rapidly changing electron 
density due to interband transition and ionization of elec¬ 
trons from valence bands causing strong spatio-temporal 
variation in excited di-electric function in non-metals. To 
couple SPPs in Cu, a series of well-characterized nano¬ 
scale grooves were machined with a micro-structured di¬ 
amond tool (overall nose radius 234 pm)using a broach¬ 
ing process (tool moves in a linear fashion to create the 
grooves repeatedly) on the 350 FG (Freeform Generator, 
Moore Nanotechnology, Inc.). This produced a nearly- 
identical family of grooves with 100 pm pitch and aver¬ 
age depths ranging from 100 pm to 100 nm, with pro¬ 
gressively narrower profile and shallower depth. For our 
experiment, the site for LIPSS formation was the out¬ 
ermost of each family of grooves, which was 1 pm wide 
and 100 nm deep, as shown in Fig. [l] The nearest neigh¬ 
bor groove was 8pm to the right (not seen in fig.). To 
avoid any influence from the other satellite grooves, the 
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FIG. 1. Depth profile (Wyko) of machined groove: 1 pm wide, 
100 nm peak-to-valley. Depth is averaged over a 10 pm length. 
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FIG. 2. Time-resolved diffraction microscopy experimental 
setup. A single pulse from the pump beam at 45° AOI forms 
LIPSS on the Cu target near a machined nano-groove. A mi¬ 
croscope objective images LIPSS formation under illumina¬ 
tion of the diffracted frequency-doubled probe beam at spec¬ 
ified time delay. The objective lens is oriented at 12° from 
normal to gather more diffracted signal. Sample is translated 
along groove structure to a fresh site for the next pump pulse. 


laser focus was positioned ~ 25 pm to the left of the nano¬ 
groove, which was over 75 microns away from the nearest 
next series of grooves on the left. 

A home-built Ti:sapphire liquid-nitrogen-cooled regen¬ 
erative amplifier produced laser pulses at 500 Hz, with 
central wavelength 773 nm, and pulse energy > 3mJ, 
which were compressed to 60 fs pulse length. Single 
pulses extracted from this train using an external Pock- 
els cell were sent into the setup shown in Fig. [2] The 
beam was split into a pump arm and a frequency-doubled 
probe arm, with the probe’s delay stage (Newmark Sys¬ 
tems NLS4-8-16-E1) variable from -0.100 ns to +1.200 ns. 
The two beams were recombined with a dichroic beam¬ 
splitter (Semrock Di02-R405) and focused onto the target 
at 45 degrees angle of incidence (AOI) in p-polarization 
using an / = 200 mm achromatic lens. The pump beam 
focus had a Gaussian waist diameter of 2wq = 43 pm 
and was filtered to 4.47 ± 0.06 pJ per pulse, correspond¬ 
ing to a beam-normal peak fluence 0.61 J/cm 2 . Before 
inserting the target, pump-probe zero-delay timing was 
determined by mixing them in a 100 pm thick nonlinear 
crystal (BBO) to produce third-harmonic signal. The 


damage site is illuminated by probe light incident at a 
known time delay, which is observed by an in-situ 20x 
infinity-corrected microscope objective (Mitutoyo) and 
imaged onto a triggered CCD camera (Basler scA1600- 
14gm). Due to their periodic nature, light scattering 
from LIPSS tends to resemble diffraction from a grating, 
while the rest of the substrate almost exclusively reflects 
in the specular direction. Then by positioning the ob¬ 
jective near 12° from target-normal, we simultaneously 
rejected the 0-order reflection and collected diffractive 
orders 2-4, resulting in high signal-to-noise for the imag¬ 
ing of the LIPSS. While plasma self-emission is often a 
nuisance in pump-probe imaging experiments, by choos¬ 
ing a low enough pump fluence in combination with a 
bandpass filter (Semrock 390/40), the pump contamina¬ 
tion is observed to be below the noise threshold of the 
CCD. 

Fig. [3] displays the data collected by this time-resolved 
diffraction microscopy method. Each false-color image 
was collected by initiating damage at a new site along 
the fabricated groove, with the probe set to the indicated 
delay time. The Gaussian maximum of the pump pulse 
can be seen by the elliptical feature, which is off-center 
from the groove by 19 pm. The corresponding plots are 
lineouts of the pixel values along the central axis, and the 
red shading indicates the groove location. LIPSS are be¬ 
ginning to form at 50-75 ps, consistent with RefP^, and 
are well developed by 200-400 ps. At longer delays, the 
LIPSS become more distorted, and by 1.200 ns the struc¬ 
ture appears very similar to its final state (not shown). 
As far as we are aware, this is the first reported observa¬ 
tion of LIPSS formation in which the formation of each 
period is resolved separately. 

We used an interferometric microscope (Wyko NT9100 
by Veeco) to profile several damage sites, and an ex¬ 
ample trace is shown in Fig. [4] The profile shows an 
elliptical central spot which is reasonably flat, in addi¬ 
tion to the LIPSS which were seen in the TRDM images. 
This flatness would explain why the central region has 
low diffracted signal in TRDM. All of these features are 
raised above the initial sample elevation, which suggests 
that at this low fluence there is minimal ablation. In¬ 
stead, we interpret that the copper is melted and amor- 
phized in these regions, and due to an expected approx¬ 
imate 4% volume increase for amorphous copper^ 3 . A 
30 nm rise then implies about 700 nm melt depth, and 
with an elliptical melt area of about 700 pm 2 we can es¬ 
timate required energy to simultaneously melt the whole 
volume Umeit ~ 0.9 pJ, which is only ~ 1/5 of the in¬ 
cident energy. LIPSS period is measured to be 2.3 pm, 
which is consistent with the prevalent SPP model for this 
AOI®. The LIPSS period is given by A = A/(A/A s ±sin6 > ) 
with AOI 0 , SPP wavelength A s = Ay/ (e' + e^)/e'e^, and 
+ or — determining backward- or forward-propagating 
SPPs respectively. In this case, = 1 for air and 
e = e' + ie" = —22.7 + il.89 the dielectric function for 
copper, which gives A = 2.45^o!i8l im due f° uncertain 
6 = 45 ± 2° from alignment. 
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FIG. 3. Time-resolved diffraction microscopy images and 
corresponding lineouts. Images are in false color with pre¬ 
fabricated groove edge marked (red line to guide the eye), and 
averaged horizontal lineout location shown (gray rectangle). 
LIPSS are beginning to form near 50-75 ps, are well-developed 
by 200-400 ps, after which they are deformed by the central 
spot evolution. Horizontal axis is labeled approximate due 
to magnification uncertainty induced by oblique microscope 
objective. 


The laser-induced damage process is generally under¬ 
stood as a fluence distribution which excites electrons in 
the interaction regime during the laser pulse, and then 
later transfers energy to the lattice via electron-ion colli¬ 
sions (i.e. the two-temperature model 14 ). After the ions 
are heated, they can propagate heat via ion-ion colli¬ 
sions and/or undergo phase transitions which change the 
surface morphology. These transport processes have not 
yet been well characterized, but determining the cause of 



FIG. 4. Study of damage site resulting from femtosecond 
pulse interaction with grooved single crystal Cu substrate, 
profiled with interferometric microscope, (a) 3D profile ren¬ 
dered with depth dimension exaggerated by ~300x. Incident 
laser k (45°) and E (p-pol.) shown. Damaged area is raised 
(crater formation absent), and mostly flat, which is consis¬ 
tent with the dark central spots observed forming in TRDM 
images, (b) 2D false color depth profile, with faint gray bar 
indicating lineout shown in Fig. [5] LIPSS period 2.3 pm. 
Peak at far right (red arrow) marks the location of the initial 
groove edge, but it clearly has also been modified by the laser. 


our time-ordered LIPSS formation result may give insight 
into the excited metal dielectric function e^t, the Drude 
collision frequency T, and the electron temperature T e . 

Our observations show a surprising result, that forma¬ 
tion of LIPSS from a pre-existing groove is time-ordered 
such that the structures form sequentially away from the 
source groove, even though the laser intensity there is the 
lowest. To explain the time-ordering, we postulate that 
the LIPSS ridge nearest to the groove had the highest 
T e , causing it to undergo ultrafast melting and expan¬ 
sion faster than subsequent LIPSS ridges. This distri¬ 
bution assumes T e gradient is opposite the to the laser’s 
Gaussian intensity gradient, due to the SPP electric field 
enhancement. It also assumes that energy stored in 
the combined E-field excites electrons, which thermalize 
within themselves via electron-electron collisions, which, 
in turn transfer heat to local ions ’adiabatically’. Based 
on these assumptions, we construct a 1-D model where 
the near surface electron energy distribution in a direc¬ 
tion perpendicular to the nano-groove is calculated by 
spatio-temporally integrating the square of the laser and 
the SPP E-field originating from the nano-groove (energy 
stored) over the laser pulse: 

/ OO POO 

I(x)dt OC / (ELaser + riEsPp) 2 dt (1) 

-oo J — oo 


where 77 is the relative field amplitude between the SPP 
and laser, and the normalized field terms are given by 
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FIG. 5. Results of Eq. [l] the ID electron heating model 
(lines), and experimental depth profile (dots). The laser heat¬ 
ing (red) in the absence of an SPP (77 = 0) is consistent 
with the central spot which forms with or without the sur¬ 
face scratch, and a relative melt threshold can be estimated 
(dashed). For SPP parameters determined by e co id^r] — 1 
(cyan), neither the hot-groove temperature gradient, nor the 
laser spot profile are reproduced. The experimental results are 
better represented using ehot (blue), with temperature gradi¬ 
ent toward the scratch, and the correct number of LIPSS are 
predicted before they are washed out by the central melt zone. 

where vp = c/sin# and v s = c/ y/ (e' + e^/e'ed ~ 0.98c 
are respectively the laser and SPP phase velocities along 
the surface, r the pulsewidth, d the laser off-center dis¬ 
tance, and wo the beam waist radius. For both fields, the 
first term is the propagating temporal profile, the second 
term is the laser spatial profile which is simply an am¬ 
plitude for the SPP, the third is the phase term, and the 
fourth term in Eq. [3] is the SPP decay due to propaga¬ 
tion. Eq. [I] can be integrated analytically (not shown), 
and the results are shown in Fig. [5] for our experimental 
parameters and various 77 and e. 

The field enhancement 77 has been calculated previ¬ 
ously for a plane wave 15 : 

2 = ESPP _ rpel _ 2 I e '\ COS#(l— R) 

' EL~r - 1) 1/2 1 ’ 

In our case this could give rise to values as high as 
77 « 9(1 — R) 1 ' 2 for room temperature copper, or 77 « 
2(1 — R) 1 / 2 if heating were to increase T by an order 
of magnitude so that Chot ~ —12 + 7 12, using simple 
equation for the dielectric function of a metal: e meta / = 
1 — uj 2 /uj{uj + iT) . Fig [3] shows that the 1 -D model quali¬ 
tatively explains our experimental observations best with 
thoti whose primary effect is reduction of the propagation 
length of the SPP by an order of magnitude, while only 
varying A by 1 ~ 2% Without this reduction of SPP 
propagation length (e co ^), the decay of the oscillations 
is dominated by the decreasing temporal overlap between 
the SPP and laser, arising due to velocity mismatch at 
this AOI. 

To arrive at this e^ot value, we first calculate T from 
the known e at room temperature, and then add in the 


electron-electron collisions T ee as the electron tempera¬ 
ture T e sharply increases due to the laser. At room tem¬ 
perature, electron-phonon collisions Y ep dominate, but 
since the ion temperature remains bas ically unchanged 
during an ultrashort laser pulse, models 1617 predict only 
a weak dependence of Y ep on T e . Then from e me t a z and 
the 300 K values we obtain 1 /Y ep = 5.14 fs. Electron- 
electron collisions Y ee are given by 17 

Tee « for k B T e < E F (5) 

with Planck constant /i, Fermi energy Ep = 7.0 eV for 
copper, and Boltzmann constant kp- While negligible 
at 300 K (1/T ee = 6.9 ps), T ee begins to dominate at 
temperatures on the scale of the electron energy in laser 
field, U p = e 2 E‘j jaser /(Am e u 2 ) = 0.53eV. Taking into 
account the off-center distance and contribution of the 
SPP field gives an effective U p - e ff = leV, which is 
still lower than the T e distribution of the blue curve 
shown in Fig. [5| which corresponds to T e = 3eV and 
T = T ee + Y ep = 1/0.5fs and best matches final com¬ 
bined LIPSS and laser melt zone profile. As ponderomo- 
tive scaling allows for electron KE up to 10 con¬ 

tributions from inter-band transitions ( 2.2 eW 9 ) due to 
conduction-band electrons collisionally ionizing valence 
electrons, could significantly change the electron density, 
resulting in higher T e , which, in turn would raise T^ on 
beyond melt threshold 14 . 

In this letter we have presented the results of time- 
resolved diffraction microscopy experiment, which for the 
first time resolved the temporal dependence of LIPSS for¬ 
mation by a femtosecond laser pulse interacting with a 
surface groove. The LIPSS were observed to form sequen¬ 
tially in time with increasing distance from the groove. 
While many aspects of LIPSS dynamics need to be stud¬ 
ied further, such as the spatiotemporal-dependence of e, 
details of the transport process, our ID model appears 
to capture most of the relevant features: LIPSS period, 
temperature gradient of LIPSS toward the groove, cen¬ 
tral melt size, and even an order of magnitude estimate 
of T e . This technique can be used widely across materials 
to understand LIPSS formation dynamics in general, and 
evolution of di-electric response of materials in particular. 
Overall, our model of electron heating due to local field 
energy density provides a reasonable explanation to the 
experimental observations. This work also motivates the 
need of further theoretical and computational efforts to 
calculate dielectric function of solids under intense fields. 

This work was supported by the Air Force Office of 
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